UVA-93-41 
December 1993 



SINGULARITY THEORY FOR W- ALGEBRA 

POTENTIALS Q 

' Jose Gaite 



Instituut voor Theoretische Fysica, University of Amsterdam, 
Valckenierstraat 65, 1018 XE Amsterdam, The Netherlands 

Q 

and 

m 

Instituto de Matemdticas y Fisica Fundamental, C.S.I.C., 
Serrano 123, 28006 MADRID, Spam. 

> 

in 
o 



On 



Abstract 



The Landau potentials of Vl^-algebra models are analyzed with algebraic- 
geometric methods. The number of ground states and the number of independent 
perturbations of every potential coincide and can be computed. This number 
agrees with the structure of ground states obtained in a previous paper, namely, 
O as the phase structure of the IRF models of Jimbo et al. The singularities 

associated to these potentials are identified. 

> 

^ 1 Introduction 

After the remarkable work of Zamolodchikov (see also Ludwig-Cardy) I] showing how 
to construct Landau- Ginzburg potentials for the minimal models of 2d CFT, there have 
been many interesting extensions and generalizations. In particular, we saw in a pre- 
vious paper that Zamolodchikov's method with some qualifications works perfectly for 
the natural generalization of the Virasoro minimal models to higher discrete symmetry, 
namely, the minimal models of W-algebras. We found the corresponding lagrangians 
and analyzed the phase structure that they entail. The chief result was the construction 
of a perturbation that produces the desired state diagrams, namely, dominant-weight 
diagrams of the corresponding Lie algebra, thus reproducing the phase structure of the 
IRF models of Jimbo et al [|| . In the case of W($) , they are triangular diagrams with 
a particular structure || . 

* Research supported in part by CICYT 



1 



However, in absence of further analysis we cannot assert that the state diagram 
mentioned above is the maximum possible unfolding of ground states. Neither do we 
know that the algebra of perturbations corresponds to the relevant fields of a definite 
model. In fact, for all we know, it might even not be finite. In comparison, these 
problems do not arise for the potentials of the minimal Virasoro models, since they 
coincide with the ADE singularities, the resolution of which is standard in the literature 

This paper is devoted to establish the unfolding of extremal points and dimension 
of perturbation algebras of W^-potentials and to prove therefore that they are well- 
defined and just right to describe the desired phase structure. We should point out that 
both features, namely, number of extremal points and dimension of the perturbation 
algebra, are intimately related and actually coincide, as will be explained below. We 
will work out in detail the low-p cases, explicitly exhibiting the perturbation algebras. 
Finally, we find that these potentials correspond to a simple series of two-variable 
singularities, of which the first two are already known. 



2 Ground states and algebra of perturbations. Gen- 
eral methods 

Let us recall the Landau potential for W? 3 \ 

V ~ a p - /L a p ~ 1 + (aa) p - 2 + c.c. 

= (aaf- 4 (a 3 + a 3 ) + (aaf- 2 . (1) 

From now on, we consider the complexified potential, letting a and a take arbitrary 
complex values and we rename them as x and y, respectively!]. Since numerical coeffi- 
cients were disregarded in obtaining ([]]) and they not play any role in this paper either, 
we will choose them as to simplify the calculations. Thus we write the potential as 

V = ^^--{xyf- A {x 3 + y 3 ). (2) 
p — 2 v ' 

It gives two algebraic equations for the equilibrium (extremal) points, 

V x = d x V = y(xy)- 5+p ({-p + 1) x 3 + x 2 y 2 + (-p + 4) y 3 ) = 0, (3) 
V y = d y V = x (x y)' 5+p ((-p + 4) x 3 + x 2 y 2 + (-p+l) y 3 ) =0. (4) 

Of course, we are only interested in the real solutions, such that y = x and the equations 
are conjugate of one another. Since they are of degree 2p — 5, one should expect either a 

^^Do not confuse with the usual notation in which x and y are real variables defined by a = x + iy, 
as in j|. 
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maximum of (2p — 5) 2 solutions or an infinite number, according to Bezout's theorem. 
In geometrical language, if the two algebraic curves intersect at a number of points 
larger than the product of their degrees, they must have a branch in common. This 
latter possibility must be excluded; it implies that the potential is not determinate, 
namely, there are higher degree terms that cannot be removed by diffeomorphisms || . 
However, it is clearly the general case in [|), where the common factor xy appears 
in both equations when p > 6. This problem was already anticipated in ||, as well as 
its solution, namely, to consider further terms in the potential. We will come back to 
it in section 5, where we construct the remaining terms. 

On the other hand, the number (2p — 5) 2 of solutions is still larger than the number 
of extremal points in the state diagram we seek to reproduce, as will be seen more 
concretely below. Nevertheless, it can still happen that the structure of those equations 
is such that the maximum number of solutions is never reached. To understand why, 
it is convenient to resort to the geometrical picture: Although the two algebraic curves 
must have the (2p — 5) 2 intersection points (counting multiplicities) in projective space, 
it can occur that some of them are permanently placed at infinity. Let us see that it 
is the case. If we homogenize the equations by adding a third variable z and we 
make z = to find the points at infinity, we inmediatly see that some of them are 
common to both curves. Furthermore, this property is not altered when we introduce 
perturbations. 

The method we will use to establish the actual number of extremal points is half 
way to solve Eqs. £|). The way one should attempt to solve a pair of equations with 
two unknown variables is by first eliminating one variable and, therefore, one equation, 
thus reducing the problem to an ordinary equation in one variable. The systematic way 
to perform this operation utilizes an algebraic object called resultant. The resultant of 
two polynomials in one variable is defined as the simplest algebraic object that vanishes 
whenever they have a common root; namely, the product of differences of roots of each 
polynomial. Since it is symmetric in each set of roots, it can be expressed in terms of 
the coefficients of the polynomials. This expression is obtained in algebra textbooks 
but will not be needed. Now, if we regard our two polynomials ([5], |^) as polynomials in 
one variable, y say, with coefficients that are polynomials in the other, x, their resultant 
is also a polynomial in x. Moreover, for those values of x that make it null the Eqs. (|3|, 
£D will be satisfied. In fact, according to our purpose, we shall not attempt to solve 
this equation in x but rather to determine the number of solutions; this can be done 
just by inspection. 

We consider now the algebra of perturbations of a potential V. This algebra is 
given by the quotient of all possible deformations by those that are induced by diffeo- 
morphisms u, 

= P[x,y] 

(v»v v y 
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Alternatively, we can consider, in place of the potential, the curves (|3]||). Each curve 
has naturally associated to it the ideal of polynomial functions that vanish on it, namely 
(V x ) or (Vy). The algebras of polynomial functions on the curves are Q x = P[x,y]/(V X ) 
and Q y = P[x,y)/(V y ), respectively. The algebra of polynomials on the intersection is 
Q. Since we are assuming that this intersection is just a finite set of points, Q is finite 
and contains as many elements as there are intersection points (counting multiplicities) 

!■ 

In conclusion, we have a suitable method to obtain the number of extremal points 
or elements of the perturbation algebra. Nevertheless, to actually determine the ones 
or the others, one is to undertake detailed computations: Respectively, to completely 
solve the equations (|3|, f|) to find the extremal points or to use these equations to 
eliminate all perturbations that they generate and see which ones are left. 



3 The Wh\ potential and its relation to the D4 sin- 
gularity 

The potential for the first model in the series, W{?\ , is well known. It was described in 
in the context of symmetric catastrophe theory. There were found the six relevant 
perturbations; with the addition of the identity, a total number of seven independent 
perturbations. Correspondingly, it was shown that the unfolding consists of seven 
points in a triangular array that constitutes the first instance of the series later ana- 
lyzed in ||. We now intend to reobtain these results within our present notation and 
objectives. 



The potential 
produces the equations 



y = x 2 y 2 _ x 3 + y 3 ^ 



V x = -x 2 + xy 2 = (6) 
V y = x 2 y-y 2 = (7) 



Given their simplicity, we can actually solve them. We thus have (0,0) with multiplicity 
4 and 3 other, (x = cubic root of 1, y — x 2 ), seven solutions altogether. It is indeed 
fewer than the expected number for two third-degree equations, which is nine. A more 
convenient way to realize this fact is to compute the resultant^, 

Resultant [V x , V y , y) = x A — x 7 . (8) 

From the geometrical point of view, the cause is clear: The curves ([], [7]) comprise two 
parabolas and two intersection points are placed at infinity. 

2 We use the Mathematica object Resultant. 
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In this case, it is also straightforward to calculate the perturbation algebra Q. It 
is clear that it contains all monomials up to the second degree, l,x,y,x 2 ,y 2 ,xy. At 
the third degree, we have to exclude right away two monomials, xy 2 ,x 2 y, because of 
Eq. (H [7]). Furthermore, we can easily obtain x 3 = y 3 , thus only leaving x 3 + y 3 . Now, 
it is not difficult to see that every remaining monomial can be expressed in terms of 
the ones just selected. For example, x 4 = xy 3 = x 2 y = y 2 . 

We must remark that, strictly speaking, x 3 +y 3 cannot be considered a perturbation 
for if it is removed from the potential, the remaining quartic term is the degenerate 
double-cusp catastrophe, an indeterminate potential. On the other hand, the quartic 
term itself can be safely removed: the cubic term constitutes by itself a well defined 
and known potential, the D4 singularity of the ADE classification. We can further 
appreciate this point by observing that the unperturbed potential has a quadruple ex- 
tremal point at the origin and three distinct others. The latter are obviously irrelevant 
for the singularity, although not for the potential as a whole, since they are precisely 
the ground states. 



4 The Wjh potential. Arnold's N\q singularity 

The next case, W^, has the potential 



3 „,3 



V -•<•// (.<""• • //•) (9) 

which has not been thorouhgly analyzed before. We know, however, that a definite 
symmetric perturbation will produce the desired triangular array of 19 extremal points, 
of which 6 are minima [[J. Proceeding as above, we consider the equations 

V x = -Ax 3 y + x 2 y 3 -y 4 = 0, (10) 
V y = -x 4 + x 3 y 2 -Axy 3 = 0. (11) 

Despite their unfriendly aspect, they are solvable. However, we must not bother to 
actually solve them, since we can get similar information from the resultant, 

Resultant [14, V y , y] = 3375 a; 16 - 27 x 19 (12) 

We see that we have indeed a maximum of 19 solutions, 16 of which are degenerate at 
(0,0), while the other 3 are again cubic roots. 

We are now to calculate the perturbation algebra Q. Right away, we can count 
in the 15 monomials up to the fourth degree. At the fifth degree, we must exclude 
two, according to [TT| ). Furthermore, less obviously, x 4 y = y 4 x. Thus we have left 



x ,y ,x y + y x. At sixth degree, there is an independent one yet, x + y , whereas 
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the remaining and any of higher degree are reducible to those already mentioned. The 
total number is 19, of course. 

We have to face now the problem analogous to that of the previous potential; 
namely, we cannot consider x y+y 4 x as a perturbation, for the same reason. Moreover, 
we have three other terms, x 5 ,y 5 , and x 6 +y 6 , that do not fit the expected perturbations 
for this potential. We are thus well advised to resort to singularity theory to further 
analyze this potential. In analogy with the previous case, we may expect the term 
x 4 y + y 4 x to furnish the appropriate singularity. Let us see whether it has the correct 
properties. First of all, we know from (|12|) that its critical point must be 16-fold 
degenerate. It is easy to convince oneself that removing the 6th degree term from the 
potential implies reducing the resultant to just x 16 . Besides, the candidate singularity 
has already been identified and named iV 16 ; it has codimension 12 and 3 modular 
parameters, associated to the monomials x 2 y 3 ,y 2 x 3 and x 3 y 3 @]. 

Therefore, we now realize that the term x 3 y 3 in the potential is just a modular 
deformation that needs not be present^. Then, the equations fllC| fLlp become homo- 
geneous of 4th degree, and indicate that two monomials of this degree do not belong 
to the perturbation algebra; let them be x 4 and y 4 . The remaining 13 monomials up 
to degree 4th account for the codimension 12. Besides, there are two other monomials 
in the perturbation algebra, x 2 y 3 = x 5 /4 and y 2 x 3 = y 5 /4:, corresponding to modular 
deformations. 

It is pertinent to point out that there is another form of the Niq singularity more 
suited to computations, namely x 5 + y 5 . However, it does not have D% symmetry (it 
has _D 5 instead) and is thus not convenient for our purposes. 



5 Wtty and the general W%\ case. The series W of 
singularities 

The potential of WfL , 



is incomplete, as was mentioned in section 2. It is easy to see that essentially we have 
only two possible additions with the necessary symmetry, x 6 + y 6 or x y (x e + y 6 ). The 
second one yields a resultant of degree 49th, showing that the number of extremal 
points is the maximum allowed by the degree of the potential. The first one yields 



Resultant[K, V^,;?/] = -1024 x 25 + 23168 x 28 - 25120 x 31 + 8944 x 34 - 1024 x 37 , (14) 



3 We must remark again though that it is essential for the potential to have absolute minima and, 
in particular, the minima that we expect. 




(13) 



6 



with 37 extrema. They account for the 10 minima, 18 saddle points and 9 maxima of 
the corresponding triangular array, as was depicted in Fig. 2}. With some effort 
one could equally obtain the 37 elements of Q but it is not worthwile. 

Let us now consider the type of singularity of this potential. The multiplicity at 
the origin is 25, according to (Hj). This is precisely the multiplicity of V = x 6 + y 6 



as can be seen at once from the equations V x = and V y = 5 As far as we know, 
this singularity has never been studied in the literature. It has codimension 18 and 6 
modular parameters. One of these modular parameters corresponds to x 4 y A , of course. 

We should attempt to generalize the previous results to arbitrary p. The first 
problem is to find out the correct way to complete the potential. To this end, we can 
appeal to the results in ||. Let us recall that we obtained there the algebraic curve 
that considered as a level curve of some potential gives the necessary extremal-point 
structure. That curve was the product 

(Z+l)/2 

I] (Io +h + Wilo) = I l +1 + Io^h + • • • = 0, (15) 
i=i 

for I odd, or 

1/2 

n(/o + h + Wih){h ~c) = I l +1 + I^h + ■ ■ ■ = 0, (16) 
i=i 

for / even. The point was that the two first terms coincide with those of the potential (|2|) 
if I = p — 3. However, it must be noticed that there are further terms that should 
be regarded as belonging to the unperturbed potential, since they remain when all 
perturbation parameters, Wi and c, are made null. Nevertheless, there is an excessive 
number of monomials in x and y; e.g., the lowest degree term in fll5|) , J^ l+1 ^ 2 } or 
in ([T5|), IqI 1 / 2 , produce when expanded a string of monomials with only two essential 
ones, x 3 ( /+1 )/ 2 -|_y 3 ('+ 1 )/ 2 or (xy) (V 3 '/ 2 + y 3l ^ 2 \ respectively. Keeping only the essential 
terms, we are thus led to propose a potential 

V = x p - 2 y p - 2 -x p - 4 y p - 4 (x 3 + y 3 )+x p - e y p - 6 (x 6 + y 6 ) + - ■ - + (x 3E ^ + y 3E ^) (17) 
for p even or 

V = x p - 2 y p ~ 2 - x p ~ 4 y p ~ 4 (x 3 + y 3 ) + x p ~ & y p ~ 6 (x 6 + y 6 ) + ■ ■ ■ + x y (x 3 + y 3 £ f 1 ) 



for p odd, which nicely generalizes that of ( |T3l) to larger p. 

It remains to be proved that this potential ( |IT| or [18]) has the properties that we 
found for those up to p = 6; namely, the maximum number of extrema is as the 

4 More generally, the multiplicity of x n + y n is (n — l) 2 . 
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corresponding triangular array and the independent perturbations can be obtained 
accordingly. Unfortunately, the form of (|l7j or [18]) as a sum of a number of terms 
proportional to p makes it unsuitable for the algebraic computation of the resultant. 
Nevertheless, we have performed the calculations up to p = 11 and checked that the 
resultant has the expected degree. The results for p = 7 to 9 are shown in the appendix. 
We believe that these partial results, beside their own interest, provide some evidence 
for the general case. 



An important result can be directly read off from (17, 18), namely, the series of 
singularities associated to our potentials: 



x 



3 ( m)/2 + y 3(m)/ 2) z = 1j3j5j 



or 

(xy) (x^ + y^), 1 = 2,4,6,-- -. 

We think that W l is a convenient name for this series^]; then W l = and W 2 = Niq. 
The latter form above, when / is even, is equivalent to ^x( 3Z / 2 )+2 _|_ ^(3//2)+2^ _ Hence, 
we see that both belong to a simple series of two- variable singularities, x n + y n . They 
have multiplicity /i = (n — l) 2 , moduli m = (n — 3)(n — 2)/2 and codimension c = 
(n + 3) (n — 2)/2 such that they fulfill // = c + m + 1 || . Among these large moduli are 
the other terms in the complete potential (|I7j, [T8D , as well as those inessential terms 
in (|T5| , jiTf ) that were discarded. 



6 Conclusions 

The Landau potentials of W3 models are analyzed in order to ascertain whether they 
fit the phase structure of the IRF models of Jimbo et al. Their sufficiency to this 
purpose has been established in a previous paper || by constructing a perturbation 
that produces the required configuration of extrema. This configuration is now shown 
to be the maximum possible unfolding. The method used only allows to consider a 
finite number of cases. The reason is that it is based on an algebraic study of the 
equations for extrema, which final step is to obtain the resultant with Mathematica. 
This is done sequentially for increasing p. It has been carried out up to p = 11 but 
only the cases up to p = 9 are displayed. The potentials for p = 4 and 5 deserve a 
more detailed study; their algebras of perturbations are also obtained. 

A necessary step to attain these results is to complete the potentials, which are 
otherwise indeterminate. This is done in general, yielding a sum of terms of decreasing 

5 We follow Arnold in giving names to series of singularities [Q. In fact, he has already used most 
letters to this effect, including W. Nevertheless, we think that it is the adequate letter here and that 
no confusion should arise since he places the other letter(s) as subscript whereas we have placed I as 
superscript. 
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degree, the last of which depends on the parity of p ( |TT| , p~8| ) . Precisely these last 
terms give the singularities of these potentials, while the others can be regarded as 
modular terms. Of course, this is the mathematical viewpoint, which disregard the 
physical significance of the various terms. An analysis along a more physical line 
should require us to resort to the actual fields associated to the monomials according 
to the identifications realized in || and is beyond the scope of the present paper. 
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A Appendix 

We show here that the resultants for the p = 7, 8 and 9 potentials precisely account 
for the number of extrema in the corresponding diagrams. Let us first calculate the 
number of extrema in the pth triangular diagram, which can be done in full generality 
by clever use of triangular numbers. This diagram has p — 2 minima on each side. 
Hence, the total number of minima is (p — 2)(p — l)/2. There is exactly one saddle 
point between each pair of minima, so that saddle points and minima form together 
a triangular array of side 2p — 5. Therefore, their total number is (2p — 5)(2p — 4)/2. 
The distribution of maxima is a little more involved; however, they can be split in 
two triangles of sides p — 3 and p — 4, the total number being (p — 3) 2 . Adding both 
numbers, we obtain 3p(p — 5) + 19 extrema altogether. In particular, 61, 91 and 127 
for p = 7, 8 and 9, respectively. 

Now we show the results of a computation of the resultants with Mathematica. The 
potential, equations and resultant are written for each case separately: 



• p = 7 




(19) 












(20) 



(21) 



Resultant [V x , V y , y] 
249221422080 



55 



12230590464 4 

1594323 X 
81199756125 



,19 



246594502656 

1594323 
8303765625 



x 



,61 



,52 



(22) 



1594323 



x 



1594323 



x 



1594323 
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• p = 8 

V = x 6 y 6 -x 4 y 4 (x 3 + y 3 )+x 2 y 2 (x 6 + y 6 ) + (23) 

V x = 3x 8 + 8x 7 y 2 - 7xV + 6xV -Ax 3 y 7 + 2xy s = (24) 

V y = 2 x 8 y - 4 x 7 y 3 + 6 x 6 y 5 - 7 x 4 y 6 + 8 x 2 y 7 + 3 y s = (25) 



Resultant [14, V^y] = 43046721 x 64 + 2027978856 x 67 - 75467102058 x 70 
+2349506627136 x 73 - 26268032446551 x 76 + 172417135016664 x 79 
+42570239554944 x 82 + 11049268755456 x 85 + 577322090496 x 88 

+990677827584 x 91 (26) 

• p = 9 

V = x 7 y 7 - x 5 y 5 (x 3 + y 3 ) + x 3 y 3 (x 6 + y 6 ) + x y (x 9 + y 9 ) (27) 

V x = 10 x 9 y + 9 x 8 y 3 - 8 x 7 y 5 + 7 x 6 y 7 - 5 x 4 y 8 + 3 x 2 y 9 + y 10 = (28) 
^ = x 10 + 3 x 9 y 2 - 5 x 8 y 4 + 7 x 7 y 6 - 8 x 5 y 7 + 9 x 3 y 8 + 10 x y 9 = (29) 

Resultant^, V y , y) = -913517247483640899 x 100 - 2700744078773949615 x 103 
+ 14078078713498713606 x 106 - 54656571881643296139 x 109 
+62947451766811208022 x 112 - 62526864228499636560 x 115 
-8771826922328936340 x 118 - 1233504732780966240 x 121 
-38877206300966880 x 124 - 67659951487955904 x 127 (30) 
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